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c« ! Abstract 

i O i 

This paper deals with the optimal path placement for a manipulator based on energy con- 
sumption. It proposes a methodology to determine the optimal location of a given test path 
within the workspace of a manipulator with minimal electric energy used by the actuators 
while taking into account the geometric, kinematic and dynamic constraints. The proposed 
methodology is applied to the Orthoglide 3-axis, a three-degree-of-freedom translational par- 
allel kinematic machine (PKM), as an illustrative example. 
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Placement optimal de trajectoires de manipulateurs pour la 
; minimisation de la consommation d'energie : 

^ ! application a l'Orthoglide 3-axes 

Resume 

L'objet de cet article est le placement optimal de trajectoires en utilisant la consommation 
energetique comme critere. Ce travail propose une mthodologie pour determiner la locali- 
sation optimale de la trajectoire dans l'espace de travail d'un manipulateur en minimisant 
la consommation d'energie electrique utilisee par les moteurs, tout en tenant compte des 
contraintes geometriques, cinematiques et dynamiques. La methodologie proposee est ap- 
pliquee a l'Orthoglide 3-axes, manipulateur d'architecture parallele a trois degres de liberte 
de translation. 

Mots cles: Placement de trajectoires, Consommation energetique, Optimisation, Manipu- 
lateur paralleles. 



1 INTRODUCTION 



Optimal trajectory planning of manipulators has been a relevant area for roboticists for 
many years. Indeed, several authors have worked on trajectory planning based on different 
optimization objectives. A review of trajectory planning techniques is given in [I]. Trajectory 
planning deals with the determination of the path and velocity /acceleration profiles (or the 
time history of the robot's joints), the start and end points of the trajectory being predefined 
and fixed in the workspace. As a matter of fact, this approach is suitable for most of robotic 
applications. A path is a continuous curve in the configuration space connecting the initial 
configuration of the manipulator to its final configuration [2]. Trajectory planning usually 
aims at minimizing the travel distance [3111], travel time [2JE1E] and/or the energy consumed 
[El El EU tH] , while satisfying several geometric, kinematic and dynamic constraints. 

Another less explored aspect of trajectory planning is the placement of a given path 
within the workspace. It aims at determining the optimum location of a predefined path to 
be followed by the end-effector of the manipulator within its workspace with respect to one 
or many given objective(s) and constraint (s). This path can be the shape of a component 
to be machined, a welded profile or an artistic/decorative profile etc. In such situations, 
the trajectory planner cannot alter the shape of the path but he/she can only play with the 
location of that path within the workspace of the manipulator in order to optimize one or 
several criterion(a). Such an approach can be very interesting in many robotic applications. 
For example, in machining, a workpiece can be better located within the workspace of the 
robot to perform a given operation more efficiently with respect to the energy consumed. 

The path placement problem has not been extensively studied in the past. Nevertheless, 
some researchers proposed to solve it with respect to various optimization objectives. Several 
performance criteria for path location problems can be considered simultaneously (multiob- 
jective) or individually, such as travel time, different kinetostatic performance indices (manip- 
ulability or conditioning number), kinematic performance (velocity, acceleration), avoiding 
obstacles, reduced wear or vibration, energy consumption etc. In the following paragraphs, 
a brief survey of the work of different researchers to solve the problem of path placement 
optimization for various applications is presented. 

Nelson and Donath [12] proposed an algorithm for the optimum location of an assembly 
task in the manipulator workspace while taking the manipulability measure as the optimiza- 
tion criterion. They considered that the location of the assembly task within the workspace 
that results in the highest manipulability is a locally optimal position for performing the 
assembly. However, Aspragathos [131 EH! considered that the manipulability index and the 
dexterity, usually quantified by the condition number of the Jacobian matrix of the manipu- 
lator, can characterize the motion ability of the manipulator, but these criteria cannot depict 
the ability of a manipulator to move in a given direction. Hence they introduced a criterion 
to characterize the best velocity performance of the robot end-effector with the path location. 
They used the concept of the orientation of the manipulability ellipsoid relative to the desired 
path and used genetic algorithm to come up with an optimal solution. 

Fardanesh et al [15] proposed an approach for optimal positioning of a prescribed task 
in the workspace of a 2J?-manipulator. Optimal location of the task is considered to be the 
location that yields the minimum cycle time for the task. In another study, Feddema [16] 
formulated and solved a problem of robot base placement for minimum time joint coordinated 



motion within a work cell. The proposed algorithm considers only the kinematics and the 
maximum acceleration of each joint in order to obtain a 25% cycle time improvement for a 
typical example. 

Hemmerle [17] presented an algorithm for optimum path placement of a redundant ma- 
nipulator by defining a cost function related to robot joints motion and limits. The proposed 
approach did not consider the path as a whole but points along the entire path, hence cost 
function considers the performance only at the node points and not the path in-between the 
nodes. 

Chou and Sadler [18] developed an optimization technique for the optimum placement of 
a robotic manipulator based on the actuators torque requirements. Pamanes and Zeghloul 
[T9l [20] considered multiple kinematic indices to find the optimal placement of a manipulator 
by specifying the path with a number of points and then assigning an optimization criterion 
to each point. The objective was to find the path location in order to have optimal values of 
all the criteria assigned to the path's points. In [21], the problem of optimal placement with 
joint-limits and obstacle avoidance is addressed. Lately, a general formulation was presented 
to determine the optimal location of a path for a redundant manipulator while dealing with 
mono- and multi-objective problems [22]. The goal of this research work was to keep the 
joint variables within their limits and to minimize the magnitude of their displacements 
[231 [211 [251 [261 [271 EH [29] . 

With a general literature survey, it comes out that although several performance indices 
are introduced or considered, there is very little emphasis on the dynamic aspects reflecting 
the energy consumption by the robot actuators. The increasing number of the robotic ap- 
plications emphasize the importance of the energy saving not only to enhance efficiency but 
also to face the world energy problems. Therefore, it is pertinent to locate well the path to 
be followed by the end-effector of a robot within its workspace in order to minimize the en- 
ergy consumed by its actuators. Accordingly, the major contribution of this paper are i) the 
minimization of the energy requirements by optimum path placement and ii) use of electric 
energy consumed by the actuators instead of treating mechanical energy relations. Hence, 
we propose an approach to optimize the location of a given path within the workspace of a 
manipulator in order to minimize the electric energy consumed by its actuators. It should 
be noted that energy optimal poses can affect various performance indices such as manipula- 
bility, dexterity, stiffness etc. In the scope of this paper, we have only considered the electric 
energy consumption as an optimization criterion (objective function) in order to highlight 
its influence on motion planning. However, others optimization criteria , such as the ma- 
nipulability, dexterity, stiffness, the motors torque required, will be taken into account in 
future work. Accordingly, we will come up with multiobjective path placement optimization 
problem that may not be convex. For that matter, we will use other optimization tools such 
as Genetic Algorithms to solve it. 

The paper is organized as follows. In section II, we propose a minimum energy path 
placement optimization problem that includes basic problem formulation, electric energy 
calculations and the algorithm proposed to solve the problem. In section III, the Orthoglide 3- 
axis: a three-degree-of- freedom translational Parallel Kinematics Machine (PKM), is used as 
a case study to implement the proposed optimization methodology and results are presented 
for rectangular test paths. 



2 PATH PLACEMENT OPTIMIZATION 



The problem aims at determining the optimal location of a predefined path in order to min- 
imize the electric energy used by the actuators. The optimization problem is subject to 
geometric, kinematics and dynamics constraints. Geometric constraints include joint limits 
and the boundaries of the workspace. Kinematic constraints deal with the maximum actu- 
ator velocities whereas dynamic constraints are related to actuator wrenches. Contrary to 
the trajectories usually defined with start and end point configurations, the entire path is 
supposed to be known within the framework of this research work. The path location can be 
defined in a similar way as to define the location of a workpiece with respect to a manipulator 
reference point. 

2.1 Path Localization 

In order to formulate and describe the problem, two reference frames are defined: i) the path 
frame T v and %i) the base frame as shown in Fig. [TJ The path frame is attached 
to the given/required path at a suitable point such as geometric center of the path. As 
J- p is attached to the path, the end-effector trajectory parameters remains constant in this 
reference frame, no matter where it is located. In other words, the path is fully defined and 
constant in T v . It can also be named workpiece frame as it characterizes the position and 
the orientation of the workpiece within the manipulator workspace. The base frame Th can 
also be called global or manipulator frame. It is attached to the manipulator base and is 
used to locate a workpiece (or jF p ) with respect to the manipulator coordinate system. The 
location and orientation of T v with respect to Tb can be defined in such a way that the 
whole path lies within the workspace. The position of T v with respect to JF b is defined with 




Figure 1: (a). Path V to be followed by the end-effector P of a manipulator, JF 6 and T v being 
the base and path frames, (b). Euler angles 

the Cartesian coordinates of the origin of T v and the relative orientation of the two frames 
is characterized by means of Euler angles. However, keeping in view the constraints of the 
manipulator wrist, Euler angles are uniquely defined in the context of milling operation with 
the parameterization given in Fig.[]Th). It allows to avoid the singularity of Euler parameters. 



As a matter of fact, any trajectory defined in T p can be transformed in the base frame Tb 
by a transformation matrix. For instance, point P, of Cartesian coordinates xp p ,yp p , zp p in 
T p can be expressed in Tb as follows: 

[p],r%[p], p W 

namely, 

[ x Pb y Pb z Pb 1 ] T Tb = b T p [ x Pp y Pp z Pp 1 ]^ (2) 

Tp being the transformation matrix from T p to Tb- Let O p (xo p ,yo p , zo p ) be the origin 
of the path frame expressed in frame Tb and (0, 0, ^) the Euler angles characterizing the 
orientation of frame T p with respect to frame Tb- Accordingly, 

cos sin 9 cos ip + sin sin -0 xo p 
sin sin # cos ip — cos sin ip yo 
cos 6* cos ip zo p 
1 

(3) 

The path placement is specified with b Tp. Let x = [xo p yo p z o p # tp] T define the 
path placement within the workspace of the manipulator, in the reference frame Tb- Conse- 
quently, the components of x are the decision variables of the optimization problem at hand. 
In the context of a general machining process like milling operation, the feature to be ma- 
chined on the workpiece is defined with respect to the frame attached to the workpiece, 
namely T p . Likewise, the machining operation conditions such as machining velocity and 
acceleration are fully defined in T p . Finally, the part to be machined is defined by the de- 
signer and located in Tb whereas the machining operation conditions and robot trajectory 
planning are defined by the production engineer with respect to the the corresponding part, 
namely T p . Here, we introduce a methodology to help the production engineer well locate 
the workpiece, namely T p , within the robot base frame Tb in order to minimize the actuators 
electric energy consumption. 



cos cos 9 cos sin 9 sin ip — sin cos i/j 
sin cos 9 sin sin 9 sin ip + cos cos i[) 
— sin 9 cos 9 sin -0 











2.2 Path Placement Optimization Problem Formulation 

The goal of this research work is to help the path planner find the best location of the path 
to be followed by the robot in order to minimize the energy used by its actuators. It can be 
formulated as an optimization problem, namely, 

"For a predefined path in T p , find the optimum location and orientation of T p with respect 
to Tb, defined by the decision variables x, in order to minimize the electric energy used by 
the manipulator actuators to generate that path, while respecting the geometric, kinematic 
and dynamic constraints of the manipulator." It can also be formulated mathematically as 
follows: 

{qu <qi< qiu 
\4i\ < 4iu (i = 1 ■ ■ ■ n) (4) 



n 

min E t = y^ Ei (x) 

i=l 



x is the path placement vector corresponding to the transformation matrix ^"Tp. E t is the to- 
tal electric energy required by the n actuators whereas Ei is the total electric energy required 
by the i th actuator to follow the path. q iy q i: Ti are respectively the i th actuator displacement, 
rate and torque, qu is the lower bound and qi U (resp. <ji U and r iu ) is the upper bound of i th 
actuator displacement (resp. rate and torque). For a given path placement vector x, these 
constraints can be evaluated by means of the manipulator kinematic, velocity and dynamic 
models. 

It is noteworthy that the manipulator geometric constraints guarantee that the whole path 
lies inside the prescribed workspace. Similarly, the bounds on actuator rates (qi U ) and torques 
(t^) ensure that the manipulator will not go through any singular configuration while follow- 
ing the path. The electric energy used by the actuators is formulated in the next section. 

2.3 Objective Function: Electric Energy 

The energy used by the motors depends on their corresponding velocities and torques. As 
a matter of fact, the electric current in the motors varies with motor velocities and torques. 
Accordingly, the motor's self-inductance phenomenon appears. The current / drawn by the 
motors and the motor electromotive potential V e can be calculated as a function of the 
required torque r and the angular velocity u of the actuators, namely, 

Ve = K e u (6) 

K t being the torque sensitivity factor or motor constant expressed in [Nm/A] and K e the 
back electromotive force constant expressed in [V^rad/sec)^ 1 ]. 
The total electric power Pt is composed of [30] : 



The resistive power loss (Joule effect): 



RI 2 (7) 



the inductive power loss: 



P L = U d l ,8) 



• the power used to produce the electromotive force: 

Pem = VJ (9) 

Accordingly, the total electric power P T can be expressed as follows: 

P T = Pj + P L + P EM (10) 



R being the motor winding resistance expressed in Ohm [Q] and L the motor inductance 
coefficient expressed in Henry [H]. 



Finally, the energy E consumed by a motor can be evaluated by integrating P? over the total 
trajectory time T, namely, 



Pt being the instantaneous electric power at time t, defined in Eq. ( flOl) . 
It should be noted that Eq. (jHJ) allows us to consider the energy used by the actuators when 
they do not move but still produce a torque to keep the manipulator at a certain stationary 
configuration (with respect to that particular direction or actuator), like resisting the gravity. 

It is noteworthy that energy calculation model presented in this section is suitable for 
the brushless motors, which are generally used for PKMs. However, depending on the mo- 
tors/derives in application, energy calculation model can be developed accordingly. 

2.4 Resolution of the Path Placement Optimization Problem 

To solve the problem, a general optimization approach is proposed as illustrated in Fig. [2j 
The approach can be summarized in three constituent elements or phases: 

1. Preparation Phase: Manipulator geometric, dynamic and electric parameters along with 
the definition of the required path are used as the known input data of the optimization 
problem. The base frame T\> is defined. The path to be followed by the end-effector 
of the robot is defined in the path frame T v . The terms of the transformation matrix 
between T v and T\, are the decision variables of the optimization problem. 

2. Evaluation Phase: At this stage, the inverse kinematic model (IKM), the inverse ve- 
locity model (IVM) and the inverse dynamic model (IDM) of the manipulator are 
determined for each set of design parameters obtained from the optimization routine. 
Accordingly, the objective function and the constraints of the optimization problem are 
evaluated. 

3. Optimization Phase: The objective function and constraints evaluated at the previous 
step are analyzed by means of the optimization algorithm. Once the constraints are 
satisfied, the objective function is tested for its optimum value. The convergence criteria 
are checked. If latter are respected, the optimization algorithm stops. Otherwise, other 
iterations are run as long as the convergence criteria are not satisfied. 

Finally, the optimum path placement is obtained by means of the position and orientation 
of the origin of T v with respect to Th defined by the terms of the transformation matrix 
between T v and namely, the location vector x*=[x£, y*o z b 0* @* i J *] T - 




(11) 



3 CASE STUDY: 

APPLICATION TO THE ORHTOGLIDE 



3.1 Description of the Orthoglide 



The Orthoglide is a Delta-type PKM [31] dedicated to 3-axis rapid machining applications 
developed in IRCCyN [32]. It gathers the advantages of both serial and parallel kinematic 




Figure 2: Flowchart of the path placement optimization process 



architectures such as regular workspace, homogeneous performances, good dynamic perfor- 



mances and stiffness. The Orthoglide is composed of three identical legs, as shown in Fig. 3(a) 



Each leg is made up of a prismatic joint, two revolute joints and a parallelogram joint. Pris- 
matic joints of the legs, mounted orthogonally, are actuated which result the motion of the 
mobile platform in the Cartesian space with fixed orientation. 

The Orthoglide 3-axis geometric parameters are function of the size of the prescribed cubic 




Figure 3: A snap and workspace of the Orthoglide 



Workspace size L wor k sp ace = 0.2 m 
Point Cartesian coordinates in T b [m] 

~o b (0,0,0) 

C (-0.027,-0.027,-0.027) 

Q+ (0.73,0.73,0.73) 

Q- (-0.127,-0.127,-0.127) 

Table 1: Orthoglide workspace parameters 

Cartesian workspace, defined by the length of the cube sides, namely, L workspace [33]. The 
base frame T b is defined with the unit vector in the direction of the i th prismatic joint, 
namely, X b , Y b and Z b , the origin O b of T b being the intersecting point of e.;. Two points 
Q + and Q~ are defined in such a way that the velocity transmission factor is 1/2 and 2 at 
these two points [31]. A cube is then constructed with Q + Q~ as its diagonal. It should be 
noted that the cubic workspace center, i.e., point C, and the origin O b of the reference frame 



T b do not coincide, as shown in Fig. 3(b) In the scope of this study, L wor \ iSvace is equal to 



0.200m. Accordingly, the coordinates of points Q + , Q and C for the said workspace are 



given in Table [TJ Similarly, the prismatic actuator bounds, p min and p ma x, can be calculated 
[22] • Table [2] shows the lower and upper bounds of the prismatic joints displacements and 
their maximum allowable velocity and torque for the Orthoglide. The geometric, kinematic 
and dynamic parameters of the Orthoglide are defined in [321 [331 EH [35] . 
Electric energy Ei used by each actuator is calculated by means of Eqs. <^ to ({Til . As 
the Orthoglide 3-axis has three 3-phase Sanyo Denki synchronous servo motors (reference : 
P3050604D), Eq. fllOp is multiplied by 3 to cater for the power consumed by the each phase 
of the motor in order to calculate the electric power Pn used by each actuator, i.e., 

P Ti =3(RI 2 + LI^- + VJ) (12) 



P^min 


0.126m 


Pimax 


0.383m 


Vj 

''max 


1.00 m.s" 1 




1.274 Nm 


l rnax 



Table 2: Orthoglide actuators parameters {i = x, y, z) 



3.2 Trajectory Planning and External Forces 

In order to apply the methodology proposed for path placement optimization, a rectangular 
test path is considered. The test path is defined by the length L and the width W of the 
rectangle, as shown in Fig. H(a). Path reference frame T p is located at the geometric center 
of the rectangle. This type of path can be the example of the generation of a rectangular 
pocket like that of Fig. H](&). The position of T p in the base frame Tb is defined with the 
Cartesian coordinates of the origin of T p , O p (xo p , yo p i z o p ) and the orientation of T p with 
respect to Tb is given by Euler's angles, as depicted in Fig. W[b). For the sake of simplicity, 
only one of the three rotation angles is considered i.e, rotation about Z^-axis while XbYb and 
X p Yp planes are considered to be always parallel. Accordingly, there are four path placement 
variables, i.e., xo p , yo P , z o P an d 0, as illustrated in Fig. [H The magnitude of the end-effector 
velocity is supposed to be constant along the path. Hence, for given path dimensions, po- 
sition vector Pjf p = [xp p yp p z Pp ] T and velocity vector Vjr p = \xp p yp p z Pp ] T in the 
path frame can be evaluated as a function of time. Figure \5\ shows the position and velocity 
profiles in T p for a 0.05 x 0.10 m rectangular path and for a constant end-effector velocity 
of 1.0 m.s" 1 . Position and velocity vectors defined in T p can be expressed in Tb by means of 
the transformation matrix defined in Eq. ([3]), namely, 
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Figure 4: Rectangular test path 
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with x = [xo p Hop z o p 0] T being the decision variables vector of the optimization prob- 
lem. For a matter of simplicity and not to deal with tangent and curvature discontinuities, 
we consider that the path is composed of four independent line segments. Therefore, we do 
not pay attention to the discontinuities between the segments. 

In order to analyze the effect of external cutting/machining forces in the generation of a 
given path, a groove milling operation is considered as shown in Fig. [6], [36J. With constant 
feed rate or end-effector velocity v p of magnitude 0.66 m.s -1 , i.e, 40 m.min -1 , the following 
components of cutting forces are considered: 
i 7 /=component in the feed direction = 10 N 
F a =component along the axis of cutting tool = 25 N 
F r =component perpendicular to Ff and F a = 215 N 

3.3 Path Placement Optimization Problem 

The path placement optimization problem for the Orthoglide 3-axis is formulated in order 
to minimize the total electric energy used by the three prismatic actuators. The kinematic, 
velocity and dynamic models of the manipulator are used to evaluate the required actuator 
displacements, velocities and torques. The constraints of the optimization problem are the 
geometric, kinematic and dynamic ones. It should be noted that within the prescribed 
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Figure 5: Test trajectory for a rectangular path of size 0.05m x 0.10m 



workspace, the Orthoglide 3-axis is free of internal collisions and that there is not any limit 
on the passive joints. Therefore, the geometric constraints for the path placement problem 
are the upper and lower limits of the prismatic joints variables. The kinematic constraints 
are the maximum velocities that the prismatic actuators can produce whereas the dynamic 
constraints are the limits of the torque/force that the actuators can produce. The kinematic 
and dynamic constraints are obtained from the catalogue, as shown in Table [2j As already 
mentioned, the decision variables are the Cartesian coordinates of the origin of T p and the 
orientation angle of T v with respect to T\>. 
The optimization problem can be formulated as follows, 



min E t = K- 



subject to: 



< T 
— 1 max 



^ Pmax 



(13) 



where x = [xo P , Hop, z Op, 4>] T ■ The subscripts x, y and z are used for three prismatic actuators 
or three Cartesian directions. p m in and p ma x are respectively the minimum and maximum 
displacements of the prismatic joints as presented in Table [1] 

The optimization problem was solved by using the MATLAB fmincon function, which is a 
general constrained optimization solver using the derivative-based search algorithms. The 
optimization process was performed with different starting points and it turned out that 
MATLAB fmincon function always converges to the same solution no matter the starting 
point. Furthermore, to study the variation pattern of energy requirements at different points 
within the workspace, a workspace discretisation is carried out with respect to the path 
placement variables and the energy is calculated for each of the discrete point for a given 
path while verifying the constraints. 
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Figure 6: Cutting forces 



3.4 Results and Discussions of the Path Placement Optimization 
Problem 

The path placement process introduced in this paper is highlighted by means of the rectan- 
gular path shown in Fig. HI As a matter of fact, the path placement process is performed 
for different rectangles with constant aspect ratio of 2, i.e., L/W = 2. With the help of the 
optimization algorithm, the location of the path corresponding to minimum and maximum 
energy consumption is obtained, i.e., the best and the worst path locations with respect to 
the electric energy consumed. Figures [7] and [8] show the location of different rectangular 
paths with the minimum and maximum energy consumption in the Orthoglide 3-axis cu- 
bic workspace. The magnitude of the energy used for both best and worst cases and the 
corresponding gain of energy is given in Table [3] and is illustrated in Fig. [91 In Fig. [9J % 
saving is the percent energy saving between the best (minimum) and the worst (maximum) 
energy consumption. It can be seen from Figs. [7| and [8] that the energy consumption is a 



Rectangular path dimensions [mm] 


F m in [J] 


Ffnax [J] 


% gain 


Width (W) 


Length (L) 


20 


40 


15.26 


44.46 


65.68 


30 


60 


22.88 


61.35 


62.71 


40 


80 


30.41 


76.31 


60.15 


50 


100 


38.55 


89.80 


57.07 


60 


120 


46.83 


102.11 


54.13 


70 


140 


56.82 


113.46 


49.92 


80 


160 


65.94 


121.17 


46.89 



Table 3: Minimum and maximum energy used for a given rectangular path 



minimum when the path is located in the vicinity of the isotropic configuration with velocity 
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Figure 7: Locations of rectangular path of different sizes (Wmm x Lmm) that yield a 
minimum energy consumption 
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X b [mm] 
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Figure 8: Locations of rectangular path of different sizes (Wmm x Lmm) that yield a 
maximum energy consumption 



and force transmission factors equal to one, i.e. point with = 0° and is a maximum 
when the path is located in the vicinity of point Q~ with = 45 °. From Fig. [9] it can be no- 
ticed that the smaller the path, the higher the energy saving. This higher gain for the smaller 
path is due to the higher range of displacement of the path within the manipulator workspace. 
In order to view the energy variation trends in the workspace, a test path of size 30mm x 60mm 
is taken and the energy required for the generation of this test path is evaluated for several 
positions and orientations. Figure [10] shows the variations in the energy used with respect to 
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Figure 9: E min and E max and percentage saving as a function of the rectangular path width 
(L = 2W) 

E vs xo p zop (0 = 0°, W = 30 mm, L = 60 mm) 
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Figure 10: Energy as a function of Xq p and zq p for a 30mm x 60mm rectangular path 



xo P and f° r a constant orientation of ° and for three different values of yo p - Figure [TT] 
illustrates the isocontours of the energy required by the motors with respect to Xo P and z/o p 
for given values of z Qp and 0, namely, zo p = and = 0°. From Fig. (TDJ it is apparent 
that the energy required is not sensitive to variations in zq p - The reason being that the 
path lies in the X^-plane and the actuator displacements along Z^-axis is not significant. 
Figures [10] and [11] also show that the energy required by the motors is a minimum when the 
path is located in the neighbourhood of the isotropic configuration and is a maximum when 
the latter is located in the neighbourhood of singularities. 
Figures [TJ] and [12] show the variations in the energy used with the path orientation in 



Energy [J] variation with xq p and yo p (zo P = 0, (f> = 0°) 
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Figure 11: Energy vs Xq p and yo P for 30mm x 60mm rectangular path 
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Figure 12: Energy as a function of xo p and yo p for different orientations (zop 



different areas of the cubic workspace. It can be seen that the energy used is usually a maxi- 
mum when <j) = 45 °. However, the energy consumption is a maximum for a path orientation 
different than 45°, for some path locations. For example, the energy required at the upper 
right corner of the workspace is higher for = 30° than <fi = 45°. Figure [141 shows the 
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Figure 13: Energy as a function of Xq p and <p f° r Hop — z Op — 
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Figure 14: Comparison of 30mm x 60mm trajectory parameters for E m i n and E max locations 



comparison of trajectory parameters of 30mm x 60mm test path for minimum and maximum 
energy locations. It makes sense that when E = E max the range of actuator displacements 
is larger than when E = E m i n . Similarly actuators experience higher values of maximum 
velocities and torques when E = E max as shown in Fig. [MT b-c). which results in higher 
energy consumption for each actuator, as shown in Fig. [MT d). These results mean that the 
actuators may reach their performance limits due to an inappropriate location of the path in 
the workspace. 
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5 CONCLUSIONS 

A methodology for path placement optimization was proposed in this paper. The electric 
energy required by the actuators to follow a predefined path is considered as the optimization 
criterion. Electric energy requirement is calculated with the help of the required actuator 
torques and velocities along with motors electric parameters. To verify the feasibility of the 
solutions, actuators performance limits such as their joint limits, maximum velocities and 
torque were used as the constraints of the optimization problem. The kinematic, velocity 
and dynamic models were used to come up with the objective function and constraints. 

The proposed methodology was applied to the Orthoglide 3-axis, a three-degree-of-freedom 
translational parallel manipulator with a quasi-cubic workspace. Rectangular test paths of 
different sizes were considered as illustrative examples. These paths are similar to those used 
to realize pocketing operations. 

The use of the electric energy instead of mechanical energy as an optimization criterion is 
pertinent. Although actuator electric energy consumption depends on the mechanical energy 
requirements, the electric energy evaluation is more comprehensive than its mechanical coun- 
terpart. General approach to calculate the mechanical energy with the help of manipulator 
velocity and dynamic models, i.e., by using actuator torques and velocities, may lead to an 
under estimation of the energy requirements in the case where actuators are experiencing 
torques with zero velocities. Besides, usual mechanical energy calculations do not consider 
the resistive energy loss in the motor windings as well as the energy loss due to the variations 
in the actuator velocities. Those variations affect the current requirements and hence induce 
electromotive forces in the actuators. Accordingly, the electric energy formulation takes into 
account all these energy losses. 

The energy required to perform a given task depends on the position and the orientation 
of the task within the workspace of the manipulator. Accordingly, some electric energy can be 
saved by properly selecting the position and the orientation of the task. Indeed, a misplaced 
task can cause excessive energy consumption and can force the actuators to go over their 
performance limits. 

For the Orthoglide 3-axis, the optimum path location is found to be in the neighbourhood 
of the isotropic configuration but there is no general rule to predict the exact optimal position 



and orientation of a task particularly for a complicated three dimensional task or for an irreg- 
ular workspace. However, a detailed analysis of the energy variation within the workspace for 
a given task can lead to the optimal position/orientation of that particular task. Numerical 
optimization algorithms are useful for such a comprehensive analysis in which all the problem 
constraints and performance measures can be considered simultaneously. 

In the future work, the path placement optimization problem will be dealt as a multi- 
objective one. For example, along with energy requirements, the manipulator dexterity and 
stiffness can be considered as optimization objectives. 
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